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My procedure was this: I would count the stones by eye and write 
down the figure. Then I would divide them into two handfuls that 
I would scatter separately on the table. I would count the two 
totals, note them down, and repeat the operation. 

Borges, Blue tigers (translated by Andrew Hurley) 


1. Introduction 

In the last decades numerous interrelations were discovered between the combinatorics 
and the theory of integrable systems. Mainly, these links involve solutions, either special 
ones, such as the Painleve transcendents [5] and solitons [9], or generic ones, such as the 
tan-function of the Kadomtsev-Petviashvili hierarchy [1]. 

On the other hand, equations themselves exhibit a certain combinatorial nature, 
due to the recurrent relations which govern the higher symmetries and conservation 
laws of integrable hierarchies. This aspect was paid less attention so far, although a 
quite simple description on the language of set partitions was known for long for the 
Burgers hierarchy [12, 13]. We reproduce this combinatorial interpretation for the sake 
of completeness and as a base for further generalizations. New results obtained in the 
paper are related to the Ibragimov-Shabat and KdV hierarchies, see table 1. In these 
cases the combinatorics becomes more complicated, since the ordinary set partitions are 
replaced by special ones which are characterized by additional restrictions. Moreover, 
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Table 1. Contents of the paper. 

Hierarchy Combinatorial objects, their numbers 


potential Burgers 
Burgers 

Ibragimov-Shabat 
Korteweg-de Vries 


Set partitions, Bell polynomials Yn, Stirling 
numbers of the 2nd kind. Bell numbers 
Set partitions without distinguished singleton 
B-type partitions, polynomials (8), H-analogs of 
Stirling numbers of the 2nd kind, Dowling numbers 
Non-overlapping partitions, polynomials (10), 
number triangle (12), Bessel numbers 


this combinatorics comes in disguise: for instance, in the Burgers case we consider the 
generating function intermediately for the higher flows, but in the KdV case we have 
to consider a formal series for the logarithmic derivative of -^-function which solves the 
Riccati equation (inversion of the Miura map, see e.g. [8]). The flows and conservation 
laws of the hierarchy are related with this generating function by simple algebraic 
relations. In the Ibragimov-Shabat case a natural choice of the generating function 
is dictated by the linearization procedure. 

Although we are not interested in ‘explicit’ formulae for the coefficients of generating 
functions here, it should be mentioned that such formulae for the potential KdV flows 
actually do exist. One of them, obtained already in [8], represents the coefficient of a 
given monomial as a certain multiple integral. Another formula obtained in [16] is of 
more combinatorial nature, but it remains very complicated. Only in the case of pot- 
Burgers hierarchy the formula for the coefficients can be considered as truly an explicit 
one. 


2. Potential Burgers hierarchy 


The pot-Burgers hierarchy is obtained from the linear heat equation hierarchy 

An = V’n (1) 

by means of the change of dependent variable ijj = e". This yields 


= e-^DAA) = iD + ui)”(l) = Ynivi, ..., u„), n = 0,1, 2,... (2) 


Here and further on we denote the derivatives as follows: Vn = D^{v), D = d/dx, 
Vtn = dv/dtn- Several first equations (2) are shown in the table 2. A meaningful 
combinatorics appears just from nothing! 

It is easy to see that Yn are polynomials with integer coefficients, homogeneous 
with respect to the weight w{vj) = j. These polynomials play a fundamental role in 
combinatorics and are known under the name of (full exponential) Bell polynomials [4]. 
An equivalent definition through the exponential generating functions reads 


^ n\ ^ ^ ^ n\ 


_ ^v(x+z)—v(x) _ 


= exp 




n=0 


n=0 


vn=l 
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Table 2. The potential Burgers hierarchy (weight w{vj) = j). 


Vto = 1 

Vtj = Vi 
Vt 2 =V 2 +vf 

Vt 3 = W3 + 3viV 2 + vf 

Vti =Vi + + S-uf) + Qv'Iv2 + vf 

vta = 'I’s + {^viVa + 10 ^ 2 ^ 3 ) + (10u3n3 + 15niv|) + lQvlv2 + 


and this immediately implies the explicit formula 


K, = 


E 


n\ 


/c 1 + 2 A:2+... +r Air ='n. 


... (r!)*^’'/ci! ... kr\ 




( 3 ) 


Its combinatorial interpretation is obvious: 


— monomials correspond to partitions of the number n; 

— coefficients of monomials count partitions of the set [n] = {1 ,... ,n} into the 
subsets (or blocks) of prescribed size. 


For example, let 

us list all set partitions for n 

= 2,3,- 

4: 


n = 2 : 

V2 

vf n = 3 : 

V3 

3viV2 

vf 


2 

1 + 1 

3 

1 + 2 

1 + 1 + 1 


12 

1|2 

123 

1|23 

2 13 

3 12 

1|2|3 


t>4 

Aviv^, 

3vl 

I>v\v2 vf 

4 

1 + 3 

2 + 2 

1+1+2 l+l+l+l 

1234 

1|234 

12|34 

1|2|34 1|2|3|4 


2 134 

13 24 

1 3 24 


3 124 

14 23 

1 4 23 


4 123 


2 3 14 


2|4|13 

3|4|12 


Recall, that each set partition is considered as unordered set (with blocks as the 
elements), that is, ordering of the blocks does not matter. However, it is often useful 
to enumerate the blocks somehow. For the sake of dehniteness, we will adopt the 
enumeration corresponding to the ordering of the minimal elements in the blocks. 

We see that the combinatorics behind the hierarchy (2) is quite simple. The 
following statement is well known, see e.g. [12, 13]. 

Theorem 1. In the potential Burgers hierarchy, the coefficient of the monomial 
is egual to the number of partitions of the set of n = ki + 2^2 + ... + rkr 
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elements into ki blocks with 1 element, /c 2 blocks with 2 elements, . .., kr blocks with r 
elements. 


Proof. One proof follows intermediately from the explicit formula (3) for the coefficients. 
However, we will not always have such a formula at hand. The following reasoning 
provides a more conceptual proof. 

Let Tln,k denotes the set of all partitions of the set [n] into k blocks and H^ denotes 
the set of all partitions of [n]. Let us consider operations 

dj ■ ^n,k t n,2+l,fc) j 1, . . . , /c, M . ^ 

dehned, respectively, as appending of the element n + 1 to j-th block or adding it to the 
partition as a new singleton. This can be visualized by the following diagram: 


k 


o 



1 -• 

1 2 


©- 


-©-© 


-©- 


o 

o 

o 

o 

o 


n 


M 

dk 


(^2 

di 


Starting from the partition {0} of the set [0] = 0 and applying operations dj,M, one 
can generate, in a unique way, any partition of [n]. Indeed, the required sequence of 
operations is uniquely recovered by deleting elements in the inverse order from n to 1. 

In the theorem, a set partition tt with ki 1-blocks, ..., kr r-blocks corresponds 
to the monomial p{7r) = Vi^ .. The differentiation D{p{7r)) by the Leibnitz rule 

amounts to replacing of Vi with Uj+i for each factor in turn, taking the multiplicity into 
account. In the partition language, this means that we add the new element to each 
block in turn. As the result, we obtain the sum of monomials p{dj{7r)) for all admissible 
values of j. Multiplication of the monomial ^(vr) by Vi gives the monomial p{M{7r)). 
Thus, the polynomials 

Pn= 

TT&Un 

are related by the recurrent relation Pn+i = {D + vi){Pn) and since Pi = vi, hence 
Pn b^('yi, . . . , Vrf). n 

A less detailed statistics is obtained if we forget about sizes of blocks and consider 
just their number in a given partition. Obviously, this correspond to summing up the 
coefficients of terms of the same degree, which gives us the Bell polynomials of one 
variable 

n 

Bn{u) = Yn{u,...,u) = {udu + uY{l) = 

k=0 
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The coefficient 


of that is, the number of partitions of [n] into k blocks, is called 


the Stirling number of the second kind [15, A048993]: 


1 1 


0 

1 






1 

0 

1 

1 





2 

0 

1 

3 

1 




5 

0 

1 

7 

6 

1 



15 

0 

1 

15 

25 

10 

1 


52 

0 

1 

31 

90 

65 

15 

1 

203 

0 

1 

63 

301 

350 

140 

21 1 

877 


By definition, = 0atn>0 and |p| = #{0} = 1. The total numbers of set 

partitions with n elements, the Bell or the exponential numbers [15, AOOOllO], are given 
by the sums of the rows: 


k=0 


n\ 


= e 




n=0 


3. Burgers hierarchy 

The right hand sides of equations (2) do not contain v and this makes the substitution 
u = Vi possible. This brings to the Burgers hierarchy 

ut^ = D{Yn{u,...,Un-i)), n = l,2,... (4) 

which is homogeneous with respect to the weight w{uj) = j + 1. Several first equations 
are written down in the table 3. What is the combinatorial interpretation in this case? 
This can be easily understood by the following example, for n = 3: 


V3 

3viV2 

vf 






U2 

3uui 


Us 3 uU2 

3ui 


Ou^ 

123 

1|23 

1|2|3 


1234 1|234 

12|34 

1|2|34 

1|2|3|4 


2 13 



2 134 

13 24 

1 3 24 



3 12 



3 124 

14 23 

414103 








2|3|14 

014143 








314140 



Certainly, renaming vj Uj-i does not change the combinatorics. The differentiation 
amounts to appending the new element to all blocks in turn, however, now we do not 
add it as a new block. Therefore, the partitions under consideration are constructed as 
in theorem 1, but we do not apply the operation M at the last step. As a result, all 
partitions n„ are mapped onto some subset of partitions fl^+i, namely, those partitions 
where the element n + 1 does not appear as a singleton. We arrive to the following 
combinatorial interpretation of equations (4). 
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Table 3. Burgers hierarchy (weight w{uj) = j + 1). 


Mti = Ml 

Ut2 = U2 + 2uui 

Mt 3 = M 3 + (3mM2 + 3Ui) + Su^Ui 

Ut4 = M4 + (4mm3 + IOU1U2) + {6u^U 2 + I2uuf} + 

Uts = M 5 + {5uU4 + 15uiU3 + IOM 2 ) + + 50uUiU2 + 15u^) 

+ (10m^M2 + + bu'^ui 


Theorem 2. In the Burgers hierarchy, the coefficient of the monomial 
is egual to the number of partitions of the set with one distinguished element into ko 
blocks with 1 element, ..., kr blocks with (r + 1) element and such that the distinguished 
element does not constitute 1 -block. 

As before, one can consider more rough statistics. For instance, setting u = 1 gives 
us the total number of partitions under consideration of the set [n + 1]: 

n 

= b;(i) = n> 1 . 

k=l 

The sequence of these numbers {2-Bell numbers) starts 

1, 3, 10, 37, 151, 674, 3263, 17007, 94828, 562595, ... 

According to [15, A005493], it can be characterized also in many other ways, in 
particular, as the number of partitions of [n] with distinguished block or as the total 
number of blocks in all set partitions of [n]. These interpretations are obvious as 
well, since the distinguished blocks can be identihed with the blocks enlarged by the 
operations dj, and these operations are applied exactly as many times as there are blocks 
in all partitions. 

4. Ibragimov—Shabat hierarchy 
4 . 1 . Recurrent relations 

The table 4 displays the sequence of point changes and substitutions between equation 
'f’ts = 'f’s and the Ibragimov-Shabat equation [10] 

Ut3 = M 3 + + Quul + 3u‘^ui. 

Although this transformation looks quite harmless, it partially destroys the symmetry 
algebra: in the variables ip, it consists of equations (1) of arbitrary order, while only 
odd order equations survive in the variables u. Indeed, the change ip'^ = s brings to 
equation st„ = ... where the right hand side is a full derivative only if n is odd: 

St„ = 2lp1pn = D{2lp1pn-1 - 2lpi1pn-2 + 2lp2lpn-3 + . . . ± '0(n-l)/2)- (5) 
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Table 4. Linearization of the Ibragimov-Shabat equation. 


V'ts = 

I = S 

Si3 = D (^S2 - ^ 


Mtg = Ms + 3 m^M2 + 9mM^ + 

I = M 

Mi3 = d(^V2-^ + 3mMi + j 


t s = gi 


f V =Wi 

M;t3 = ws - + 3 m;im; 2 + w? 



q=e‘ 


.2w 


In the analogous equation for even n the term remains outside the parentheses, 
that is, St^ ^ ImD, and therefore the further substitution s = qi leads out of the class 
of evolutionary equations. The structure of odd flows is described by the following 
statement. 

Statement 3. Let us denote Dt = dt-^ + z'^dt^ + z'^dt^ + ..., A = A{z) = oq + + 

02 ^^ + A = A{—z), then the Ibragimov-Shabat hierarchy is equivalent to equations 


Dt{u) = ^D{AA) = 1(A -A)- uAA, 
z(B + u^)(A) = A — u. 

Proof. Let us consider the generating function 
M/ = -0 + + 'ljj2z'^ + . . . 

and set T = y/2e'^A. Equation (7) follows from the relations 


( 6 ) 

(7) 


zD{^) = T — -0, fj = y/ql = \/2e^'^wi = y/2e^u. 

Next, let T = T(—^), then (cf (5)) 

D{^^) = z-\'^ - - t/j) 

= z~^^p{^ - T) = + fj^z^ + ...) = 2'4)Dt{'ip) = Dt{s). 

Applying D~^ yields TT = Dt{q) = 2e‘^'^Dt{w), wherefrom 


2uDt{u) = Dt{v) = DDt{w) = = D{AA). 


Second equality in (6) follows after elimination of derivatives by use of (7). □ 

Equation (7) is equivalent to recurrent relations 

ao = M, On = an{u, . . . ,Un) = {D + U^){an-l), 77, = 1, 2, . . . (8) 


which are our object of study. Let us try to hnd a combinatorial interpretation for this 
recursion. 

Several hrst polynomials are presented in the table 5. Given these data as 

a prescribed statistics, our goal is to hgure out a dehnition of the corresponding 
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Table 5. Polynomials a„ (weight wluj) = 2j + 1). 


ag = u 

ai = ui + 

a2 = U2 + 4u'^ui + 

as = ug + {5 u^U 2 + 8uul) + 9u'^ui + vJ 

04 = M 4 + ( 6 m^U 3 + 26uuiM2 + 8ul) + (14m^U 2 + 44u^u^) + IGu^Ui + u® 

05 = M5 + (7m^U4 + 38uuiU3 + 26uu2 + 5OU4U2) + (20m^M3 + 170m^uiM2 + 
+ (30u®M2 + l^Ou^uf) + 25u®mi + 


combinatorial objects, that is, to solve a kind of inverse problem of the enumerative 
combinatorics. In contrast to the Burgers hierarchy case, here we do not know an 
explicit formula for the coefficients, but this is not too important, the main problem 
is to guess what are the objects which we are counting. An invaluable aid in such an 
ill-posed problem may be obtained by comparison with the known data collected in the 
Encyclopedia of integer sequences [15]. Let us pass to the less detailed statistics by 
gluing together terms of the same degree. Polynomials of one variable a„(M ,... ,u) = 
{udu + u‘^)"'{u) contain only odd powers of u and their coefficients constitute the triangle 


1 







1 

1 

1 






2 

1 

4 

1 





6 

1 

13 

9 

1 




24 

1 

40 

58 

16 

1 



116 

1 

121 

330 

170 

25 

1 


648 

1 

364 

1771 

1520 

395 

36 

1 

4088 

1 

1093 

9219 

12411 

5075 

791 

49 1 

28640 


which turns out to be known: according to [15, A039755] this is the triangle of analogs 
of Stirling numbers of the second kind for the so-called B type set partitions. The 
sums of numbers in rows, that is, the total sums of the coefficients 0^(1,..., 1) form 
the sequence [15, A007405] of the Dowling numbers, or B-analogs of the Bell numbers. 
This gives us a broad hint at a possible connection between polynomials (8) and B type 
partitions. This guess is proved in the next section. 

4-2. Generating operations for type B set partitions 

Special classes of set partitions appear when one takes into account some additional 
structure of the set. Set partitions of B type (or signed set partitions, Z 2 -partitions) 
[6], see also e.g. [2, 17] make use of the reflection j —j. 

Definition 1. A partition vr of the set {— n, . .. ,n} is called the Bn type partition if; 

1) TT = —TT, that is for each block (3 E ti also —(3 G vr; 

2) TT contains only one block ttq G tt such that tto = —ttq. 
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We will denote 11^ the set of all such partitions and those partitions which 
contain k block pairs. 

In a brief notation for B type partitions, the negative elements of the 0-block are 
omitted, and only that block of each pair is displayed for which the element with minimal 
absolute value is positive; the minus signs are denoted by over bars. For instance, in 
this notation the partition —5, —4|—3, 0, 3|— 2 ,1|—1, 2|4, 5 is represented as 03|12|45. A 
graphical representation is clear from the diagram 

0-0 0-0 

•-• •-o 

•-•-• —> •-• 

•-• 

0-0 

Now let us dehne the map p from 11^ into the set of monomials on the variables 
Uj. Let 1/91 denote the number of positive elements in the block /3: 

|/3| = #{i e /3 : i > 0}. 

It is clear that the number of negative elements in the block is \f3\. Let a set partition 
71 G consists of 0 -block ttq and block pairs tti, tti, ..., Tr^, tt^, such that the element 
of TTj with minimal absolute value is positive. For such a partition, let 

p(7r) = ■ UlTril-lltls-jl • • • M|jrj.|_l'U|s-j.|. 

As an example, let us write down Ilf partitions, collecting together all partitions 


corresponding to the 

same 

monomial: 

U 3 

5u‘^U2 

8 uuf 

Qu^ui 

0123 

0|123 

0|123 

0|12|3 1|2|3|4 


0 123 

0 123 

0 12|3 


012|3 

01|23 

0 13|2 


013 2 

01 23 

0 13|2 


023 1 

02 13 

0 23|1 



02 13 

0 23|1 



03 12 

01|2|3 



03 12 

02 1 3 




03 1 2 


The resulting polynomial is exactly 03 . The following theorem demonstrates that this 
is not just a coincidence and the polynomials a„ are, indeed, the Z 2 -analogs of the full 
exponential Bell polynomials W. 

Theorem 4. The polynomials (8) are equal to 

an = 

Trens 
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Proof. Let us denote the sum in the right hand side Pn- Obviously, p^ = u = Oq, 
so we only have to prove that pn satisfy the same recurrent relations as a^, that is, 

Pn = {D + vf){pn-i). 

Notice, that deleting of elements in from any Bn type set partition gives us a -Bn-i 
type set partition. Therefore, 11^ is constructed from ^n_i by adding in in all possible 
ways. It is easy to see that this is done by the following operations: 

do ; ^n-i,k insertion of both elements in into 0 -block; 

dj : Yln-i,k J = 1 ) • • •) insertion of in into blocks iTT^-; 

dj : n^_i J = 1 ) • • •) insertion of in into blocks 

^ —)■ ^n,k+li adding of the new block pair {—n}, {n}. 

Starting from the trivial partition of the set {0} and applying these operations, one can 
obtain, in a unique way, any B type set partition. Let us keep track of the monomial 
p(7r), TT G I^n-ik under these operations: 

do : the factor n| 7 ro| is replaced with n|jro|+i; 
dj : the factor is replaced with 

dj : the factor u\^.\ is replaced with 
M : two new factors u are added. 

Therefore, application of all possible operations maps the monomial p{'k) to the sum of 
monomials (D-|-n^)(p( 7 r)). □ 


5. Korteweg—de Vries hierarchy 
5.1. Recurrent relations 

Let us recall (for a proof, see e.g. [8]) a computation method of the KdV conservation 
laws and flows, based on solving of the Riccati equation 

D{f) + f = X-u, X = zyA 


(9) 


by the formal power series 


^ , fl{u) , f2{u,Ui) , , fn{u,...,Un-l) , 

Jo)— o”*” 2 +•••+ -T ■ ■ ■ 

2 Z Z^ 

Equation (9) is equivalent to the recurrent relations 

n—1 

fl = U, fn+l = D{fn) + ^ fsfn-s, U = 1,2, . . . 

s=l 

which will be the main object of our study. Several polynomials /„ are written down in 
the table 6 . The flows are computed from the polynomials with odd subscripts: let 

g2m-l 1 


( 10 ) 


/ X _ 1 91 93 

~ 2z z^ z^ 




2ifiz)-fi-z)) 
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Table 6. Polynomials /„ (weight w{uj) = j + 2). 


fi=u 
f2 = Ml 
fs =U2 + 

/4 = M3 + 4mmi 

/s = M4 + {6uU2 + 5 uf) + 2u^ 

/e = M 5 + {8uu 3 + I8U1U2) + IGu^ui 

f 7 = UQ + {IOUU 4 + 28miM3 + I 9 M 2 ) + (30u^U2 + 50uu1) + 5u^ 


which is equivalent to recurrent relations 

m 

9l = '^1 92m+l = /2m+l + 2 92s-lf2m-2s+l-, TU = 1, 2, . . . 

s=l 

then the KdV hierarchy reads 

Ut 2 m +1 = D{g 2 m+i) = U 2 m +1 + • • •, m = 0,1, 2,... 

Moreover, polynomials (10) with odd subscripts serve as common conserved densities 
for all these flows. 

One interpretation of the polynomials /„ can be seen intermediately from the 
recurrent relations (10). Let ns consider expressions builded from the variable u 
and operations M{a,b), dj{a), 1 < j < deg a where deg a is equal to the number of 
instances of u in a. Such expressions can be called ‘unexpanded monomials’. For any 
expression ip its value expand(<p) is computed according to the following rules: 

— independently on the order of operations, all dj are applied before M; 

— the action of dj{a) amounts to replacing of j-th instance of Ui in a with Uj+i {u is 
identified with uq, as usual); 

— M(a, b) is replaced by the product ab. 

Let denote the set of all expressions with the total number of symbols u, d, M equal 
to n. For instance: 



unexpanded monomials 

expanded monomials 

n = 1 

u 

u 

n = 2 

di{u) 

Ul 

n = 3 

di{di{u)), M{u,u) 

U2, V? 

n = 4 

di{di{di{u))), 

U 3 , 


di{M{u,u)), d 2 {M{u,u)) 

UUl, UUl 


M{di{u),u), M{u,di{u)) 

UUl, UUl 
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Theorem 5. The number of different expressions builded from symbols M, dj, u with the 
same monomial as their value is equal to the coefficient of this monomial in polynomials 
fn- In other words, 

/n = expand((^). (11) 

Proof. Any expression from $n+i; n > 0 is either of the form dj{a) where a G 
1 < J ^ deg a or of the from M(a, b) where a G b G Taking into account the 

obvious properties 

deg a 

expand(dj(a)) = H(expand(a)), 
i=i 

expand(M(a, 6)) = expand(a) expand(6), 

this implies that polynomials (11) satisfy the recurrent relation (10). □ 


This interpretation is fairly intuitive, but it is desirable to compare it with 
something more standard. As before, let us pass to polynomials of one variable by 
collecting together terms of the same degrees. This brings us to a number triangle 
which apparently is not in the OEIS: 


1 1 
1 1 


1 

1 




2 

1 

4 




5 

1 

11 

2 



14 

1 

26 

16 



43 

1 

57 

80 

5 


143 

1 

120 

324 

64 


509 

1 

247 

1170 

490 

14 

1922 

1 

502 

3948 

2944 

256 

7651 

1 

1013 

12776 

15403 

2730 42 

31965 


( 12 ) 


Nevertheless, the sequence of coefficients sum totals turns out to be known: /„+i[l] is 
equal to the number of non-overlapping partitions of the set [n], or the Bessel number 
Bf [15, A006789]. Notice, that identifying of all Uj results in the Riccati equation 
udu{f) -\- p = X — u which is, indeed, equivalent to the Bessel equation. Moreover, 
one can see in the triangle the Euler numbers [15, A000295], the Catalan numbers [15, 
A000108] and powers of 4. 


5.2. Generating operations for non-overlapping partitions 

This class of set partitions was introduced in [7], see also [3, 11]. Its dehnition engages 
the order relation on the partitioned set [n] = {1 ,..., n}. 


On the combinatorics of several integrable hierarchies 
Definition 2. Blocks a and /3 of a set partition vr overlap if 


13 


min a < min (3 < max a < max f3. 

The set partition is called non-overlapping (NOP) if any two blocks in it do not overlap. 
All NOPs of the set [n] will be denoted fl*. 

The interval [min a, max a] is called the support of the block a. The above dehnition 
of NOP is equivalent to the property that supports of any two blocks either do not 
intersect or lie one in another. The left diagram below shows overlapping blocks and 
the right diagram shows non-overlapping ones: 

@-@ 

®-#-# ®®® 

m -•- m -•- m 

Remark 1. A neighbour class of non-crossing partitions is characterized by a more 
restrictive condition which forbids the pattern ai < < a 2 < /32 for any elements of 

any two blocks. It is under active study in combinatorics as well, moreover, it makes 
sense to combine such types of restrictions with symmetries like the reflection for the B 
type partitions, see e.g. [14]. It is an open question, whether some integrable hierarchies 
may be associated with such kind of objects. 

Some simple properties of NOPs are the following. 

— At n = 0,1, 2, 3 we have 11* = II^ and there is only one overlapping partition 
13|24 in n4,. 

— Singletons do not overlap with any block. 

— NOPs containing only doublets can be easily identihed with the balanced sets 
of parentheses: 

0-0 

^ (()())() 

•-• 

The last property explains where from the Catalan numbers appear in the triangle 
(12). The deletion of differentiation in equation (10) brings to the recursion for the 
‘dispersionless terms’: 

n—1 

fl=u, fn+1 ='^fsfn-s M,0,M^0,2n^,0,5M^,0, ... 

S=1 

In order to establish a correspondence with the general polynomials /„, let us 
identify the variable u with the set partition {0} and define the action of the operations 
M and dj on the NOPs, in such a way that expressions <hn+i be in a one-to-one 
correspondence with 11*. 

Degree. Let deg tt = fc if tt contains k — 1 multiplets (blocks with more than one 
element). 
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Operation M. Let p G 11*, a G H*. Denote by (cr)r+i the partition of the set 
{r + 2,r + s + l} obtained from a by adding r + 1 to each element, and dehne 

a) = p U {{r + 1, r + s + 2}} U G 

This can be illustrated by the diagram 

o 


o 

0-0 


X 


o 

0 - 0—0 


» 

o 

0-0 




In particular, if p = {0} then (cr)i is bounded by the doublet {1, s + 2}, and if a = {0} 
then the doublet {r + 1, r + 2} is appended to p. Notice that deg M (p, a) = deg p deg a. 


Operation dj. It consists of adding one element n + 1 to tt G If*. If j = 1 then 
the element is added just as a singleton. For 1 < j < k = degvr, the operation 
requires a detailed description. Let us denote p 2 , ■ ■ ■, Pk all multiplets in vr, ordered by 
increase of their minimal elements. Assume that all blocks with support containing pj 
are enumerated by a sequence ji < ... < js = j ■ Let us divide each of these blocks into 
the left and right parts with respect to m = max/ij-: 

% = {* e i < m}, ={ie : i > m} 
and form the new blocks 


hii = hji U {m,n +1}, r = 2 ,..., s 

as shown on the following diagram. The rest blocks of the partition do not change under 
this operation. 


j = js 

32 

ji 


9 - 


9 - 


9 - 




■ 9 - 


■9 - 9 - 


m n + 1 



djf 


9 - 


9 - 9 - 

9 -» 


■9 - 9 - 




■ 9 - 


■ 9 - 



■9 


Theorem 6. The operations M, dj generate any non-overlapping partition, in a unigue 
way. 

Proof. The last operation bringing to a given partition is uniquely dehned by 
consideration of the block f3 containing the maximal element of the partition. If it 
is a singleton, then the last operation was dp, if it is a doublet, then it was M; if it is 
a multiplet, then the operation was dj where j is the maximal number such that the 
support of multiplet pj contains the last to the end element of (3. In each case, applying 
of inverse operation brings to NOPs with lesser numbers of elements. □ 
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Taking the theorem 5 into account, the established bijection allows to associate a 
certain monomial with each NOP, although not in a quite effective way, because we hrst 
have to build an exression G corresponding to vr G n*_^ and then to compute 
expand (</?): 

G4 n* 1 

expand f ^ 

fn 

Nevertheless, it is easy to trace at the degree of monomial under this correspondence; it 
is one more than the number of multiplets in the partition. This gives us the following 
interpretation of the number triangle (12). 

Corollary 7. The number of NOPs of n elements containing k multiplets is egual to 
the number in the n-th row and k-th column of the number triangle (12), starting their 
enumeration from 0. This number is egual to the coefficient of in the polynomial 
Fn+i{u) = fn+i{u,... ,u) defined by the recurrent relations 

n—1 

Fi = u, Fn+i = udu{Fn) + ^ FsFn_s, n = 1, 2,... 

S = 1 


6. Conclusion 

We have established a relation between several classes of set partitions and generating 
functions for integrable hierarchies. Hopefully, this observation may turn useful for 
both theories. Of course, we have too few examples at the moment to make far- 
reaching conclusions. A conjecture is that each integrable hierarchy has an underlying 
generating function which may be interpreted as statistics for some kind of combinatorial 
objects (possibly unknown). As further steps, it would be interesting to reveal the 
combinatorics associated with the niKdV equation, KdV-like equations of 5-th order, 
nonlinear Schrodinger equation and so on. 

On the other hand, the objects studied in the combinatorics are so plentiful and 
diverse that it seems doubtful that any one can be associated with an integrable 
hierarchy. In all likelihood, this is a very special property, so it would be interesting 
to understand what is the integrability intermediately in combinatorial terms (rather 
than of the level of generating functions). In particular, one can try to obtain a proof of 
commutativity of the flows of a hierarchy based on their combinatorial interpretation. 
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